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, Abstract 

0^ i We present methods for obtaining new solutions to the bispectral problem. 

^ I We achieve this by giving its abstract algebraic version suitable for generaliza- 

^ i tions. All methods are illustrated by new classes of bispectral operators. 
T— I \ 

T— I ! 
O ; 

^ ■ Introduction 
\o ' 

' The bispectral problem of J. J. Duistermaat and F. A. Griinbaum Q consists of 

bi)' finding all bispectral ordinary differential operators, i.e. operators L{x,dx) having a 

i family of eigenfunctions ip{x^z)^ which are also eigenfunctions for another differential 

qh| operator K{z,dz) in the spectral parameter: 

.> ; L(x, dx)ij{x, z) = f{z)^{x, z), (0.1) 

X ■ A{z, d,)ip{x, z) = e{x)i;{x, z). (0.2) 

■ " " ' In we constructed large families of solutions to this problem, generalizing all 

previously known results (cf. [^, 17, 13|, etc.). We obtained them as special 



Darboux transformations (called "polynomial") from the most obvious solutions - 
the (generalized) Bessel and Airy ones. We recall that the Bessel (respectively - 
Airy) operators have the form L = x~^ P{xdx) (respectively L = P{dx) — x) for 
some polynomial P of degree A^. 

For fixed tjj, all operators L{x,dx) for which ip is an eigenfunction form a com- 
mutative algebra A^, called the spectral algebra. Then Spec^^ is an algebraic curve 

- the spectral curve. The dimension of the space of eigenfunctions ip (= g.c.d. 
ordL, L G A^p) is called a rank of A^. Following G. Wilson |17] we call the spectral 



algebra A^ bispectral iff there exists an operator A{z,dz) satisfying (|0^ 

The purpose of this paper is to present general methods for constructing bispec- 
tral operators (respectively - algebras). Some of our statements are only abstract 
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versions of concrete results |^ [T^, |3|, |T4|. But adopting the algebraic point of view 
enables us to clarify the ideas and to give many classes of new examples of bispectral 
operators. Moreover, being purely algebraic, some of the methods work also for dif- 
ference or (^-difference analogs. The latter case is illustrated here by many examples 
but more complete theory in the spirit of will be considered elsewhere. 

In the process of our work |^ we substantially used the bispectral involution b 



introduced by G. Wilson |17|. It is, roughly speaking, the map exchanging the roles 



of the arguments in the wave function ^(x, z): 

bijj{x, z) = ip{z,x). (0-3) 

The main idea which we exploit here is to consider the bispectral involution as an 
anti-isomorphism of algebras of ordinary differential operators (cf. |jl^), or more 
generally - arbitrary rings. 

Our first construction of new bispectral algebras from old ones is based on ex- 
hibiting anti-isomorphisms of algebras as compositions of already known ones and 
inner automorphisms. In particular, we find new bispectral algebras of arbitrary 
non-prime rank with both spectral curves (in x and in z) isomorphic to C. 

The second construction is an abstract version ~ which again uses only the bis- 
pectral anti-isomorphism - of the bispectrality theorem from Analizing the latter 
has led us to the notion of bispectral Darboux transformations which is applicable 
to other (than differential) operators. 

At the end we consider sequences of Darboux transformations depending on 
integers parametrizing the elements of the sequence. Then allowing the parameters 
to be any complex numbers we again get new solutions to the bispectral problem. 

1 Bispectral triples 

Here we present a general formulation of the bispectral problem which works not 
only for differential but also for difference or (/-difference modifications of the original 
problem of Q. 

Let TZ and TZ' be associative algebras over a field k and let 7W be a left module 
over both of them. Let /C and /C' be fields such that k C IC C TZ and k d fC' <ZTZ' . 

Definition 1.1 We call an element L ^ TZ bispectral iff there exist ^ (z A4, A TZ' , 
f e IC', 9 e fC such that 

LV = /V', (1.1) 

AV' = eij. (1.2) 

□ 

In order to have non-trivial problem we assume that 

Li{j = implies L = for L G 7^ (or TZ'). (1.3) 



Let us fix 7/^ G satisfying (LI) and (|1.2| ). We are interested in the equation 
(cf. 0) 

= W (1-4) 
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for P £TZ,Q gTZ'. We put 



= {P en\3Q for which (1.4) is satisfied}, (1.5) 
B'^ = {Q GTZ' \3P en for which is satisfied}. (1.6) 



We shah assume that the actions of TZ and TZ' on A4 commute and that (1.3) holds. 
Then B^ and B'^ are associative algebras over k without zero divisors. Obviously 
(|1.4D defines an anti-isomorphism 



b:B^^B'^, h{P) = Q. (1.7) 
Introduce also the subalgebras 

lc^ = B^n K, /c; = n /c', (1.8) 

A^ = b-\JC'^), A'^ = b{}C^). (1.9) 

Then Ai(, is a commutative algebra isomorphic to JC'^ (the spectral algebra). With a 
tilde we shall denote the fields of quotients of these commutative rings, e.g. /C^ C 
/C^, C /C. Then b extends to isomorphisms /C^ — > A'^ and A^ — JC'^. 

Definition 1.2 We call the triple {B^,B[p, b) a bispectral triple iff both /C^ and /C^ 
contain non-zero elements. □ 

Obviously if {B^, B'^, b) is a bispectral triple then any element L G A^i, is bispectral. 

Example 1.3 (i) Let TZ = TZ' he two copies of the rings of differential operators 
of one variable with rational coefficients and let be a module of differentiable 
functions '4>{x, z) (or formal power series) on which IZ acts on the variable x and TZ' 
acts on the variable z. Put ip{x, z) = 5{x — z) (Dirac delta function). Then B^ and 
B'^ are two copies of the Weyl algebra of the differential operators with polynomial 
coefficients and b is the formal conjugation *, i.e. b{x) = z, b{dx) = —dz- 
(ii) P, 14 1 Let La be the (generalized higher) Airy operator 



N-l 



i=2 

and let ^{x) be a solution of the equation = 0. Put ip{x, z) = $(x -|- z). Then 
B^ = B[p are again two copies of the Weyl algebra but now b is defined via 

b{x) = La, b{d) = d. 

Then K,^ = C[x], A^ = C[Lq,]. 

(iii) Let Lp be the (generalized) Bessel operator 

Lp = x-''{D-I3i)---{D-(3n), pec"", 

where D = xd. Let ^{x) be a solution of the equation L^$(x) = <&(x) and ip{x, z) = 
^{xz). Then for "generic" P (see |p) B^ = B'^ is generated by x^,D,Lp; now 
b{x^) = Lp, b{D) = D, biLp) = x^ and /C^ = C[x^], A^ = C[L^]. □ 
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Example 1.4 (i) Let TZ and TZ' be two copies of the algebras of g-difference opera- 
tors with rational coefficients, i.e. operators of the type 

+ ai{x)D^-^ + --- + aN{x), 

where Dg is the dilatation operator, Dqf{x) = f{qx), g E C, and ai{x) are rational 
functions. Let ip{xz) = expq(xz) be the g-exponent (see e.g. [11]) 



expq(x) = 



CO ^ 

X^' 



where {q; q)n = {I — q){l — q^) ■ ■ ■ {I — q^), {q; q)o = 1- Let B = B' he two copies of 
the g-deformed Weyl algebra of operators 

'Yaj{x)d^g, 

where dg = {Dg - l)/x, aj{x) G C[x]. Put b{x) = -dg, b{dg) = -x. Then {B,B\h) 
is a bispectral triple. 

(ii) Let Lg^p, [3 S C^, be the g-deformed Bessel operator: 

One can find solution to the equation 

of the form 

oo 

1p{x) = X^'' ^CjX-' , Oo = 1. 
j=0 

Define the algebra B to be generated by Dg, x^ and Lf^^g and B' = B (in B' the 
variable will be denoted by z). Then the anti-izomorphism will be 

Lp,g ^ (LIO) 
b: x^ ^ Lp^g (1.11) 

Dg ^ Dg (1.12) 

with a joint eigenfunction tp[x,z) = ip{xz). □ 



Example 1.5 Let ip{x,n) = 2~"Hn{x), x € C, n G Z, where Hn are the Hermite 
polynomials (see e.g. Q) for n > and Hn = for n < 0. Put B to be spanned by 
X, dx and B' to be spanned by n, T, T~^, where T is the shift operator on functions 
defined on Z : {Tf){n) = f{n + 1). Then an anti-isomorphism b can be given by 

b: x^(r + ^T-i) (^_;L3) 
d^riT-^. (1.14) 

It is easy to check that the Hermite operator — 2xd is bispectral. One can repeat 
this construction for any set of orthogonal polynomials. □ 

This example indicates that a general theory should work in various situations. 
In the next section we will see how to build new bispectral triples from a given one. 
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2 Automorphisms of algebras of differential operators 



In this section we will assume that the algebras B and B' consist of differential 
operators. The main tool for constructing new bispectral operators will be the 
following simple observation. 

Proposition 2.1 Let {B,B',b) be a bispectral triple and let L £ B. // adL is a 

locally nilpotent operator adL : B ^ B then a = e^^ is an automorphism of B and 
{B, B', boa) is a bispectral triple. 

(Here locally nilpotent means that for each M B there exists n such that (adL)"M = 
0.) The proof of Proposition 2.1 being obvious is omitted. □ 

An important example of an operator L acting locally nilpotently on B is any 
bispectral operator Q. This can be proved by using the bispectral anti-isomorphism 
b and the fact that the function / = b{L) acts locally nilpotently on B'. 

Remark 2.2 (i) Obviously, exchanging the roles of B and B' one can use also the 
anti-isomorphism 62 = e^*^^ 06, where L' acts locally nilpotently on B'. We will use 
Proposition |2.1| in both versions. 

(ii) In the definition of a bispectral triple a decisive role is played by the function 
G Ai via which the anti-isomorphism is defined. Proposition |2.1| defines a 
new anti- isomorphism bi in terms of the old function ip (see (|1.4|)). But to have 



really a bispectral triple in the sense of Definition |1.1| we face the difficult task to 
define a new function ij) to satisfy Pijj = bi{P)ijj. The most natural choice ij) = e^tj) 
often makes no sense. In the examples below after defining the bispectral triple 
using ifj the new function (or formal power series) ip can be determined using other 
arguments. □ 

Although trivial the above Proposition is a source of new bispectral operators. 

Example 2.3 (i) Let B and B' be two copies of the Weyl algebra and let b{x) = dz, 
b{dx) = z, -0 = e^^- For any polynomial p(^) the operators adp(x) and &(lp{dx) act 
locally nilpotently on B. Then if we put ui = e'^^'^^^ and 61 = 6 o ci we have 

hi{x) = b{x) = dz, (2.1) 
bi{dx) = b{-p'{x) + dx) = z- p'idz). (2.2) 

In this way we get a new bispectral triple {B, B', bi) with a new function ip{x, z) = 

gP'(a:)gX2 

(ii) One can apply the automorphism a2 = e^'^'^^^'^^ to the triple {B, B' , bi). Then 
we get 62 = ^1 o C2 , which gives 

b2ix) = 6i(x + q'idx)) =dz + q'{z - p'{dz), (2.3) 
b2{d..) = bi{dx) = z-p'{dz). (2.4) 

For the new function ip2ix, z) we must have 

dxil>2 = {z-p'{dz))^2, (2.5) 
x^2 = {dz-c^{z-p'{dz)))ilj2. (2.6) 
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This implies 



(2.7) 



which gives a new non-trivial example of bispectral operators, provided ip2 exists. 
One can construct a solution 1^2 using Fourier transformation with respect to both 
z and X on (2^) and on 

Because /C^j — C^[^]> ^!/)2 
C. □ 



xi;2 = {dz + q'{d^))i;2. (2.8) 
'[z] we see that both spectral curves are isomorphic to 



Remark 2.4 The automorphisms considered in Example 2.2 generate the entire 
group of automorphisms of the Weyl algebra 0] . This result leads to a classification 
of all solutions to the bispectral problem such that = B'^ is the Weyl algebra. 

In general, we propose the following approach to the classification problem. 
Given solution ijj to the bispectral problem is first reduced to a one with a "better" 
B^ (perhaps by a "bispectral Darboux transformation" - see the next section). In 
all known to us cases this B^ is either the Weyl algebra or {x^ , D,Lp) (see Example 
pT^ ). Then the second step is for a given ring B to find all bispectral triples {B,B',b). 
□ 



One can easily check that the operators in the r.h.s. of ( |2.5| ) and (2.6) satisfy the 
so-called string equation [-P, Q] = 1, which plays an important role in 2d-quantum 
gravity |5|, H, Q (see also |l5|). The corresponding to the r.h.s. of ( p3| ) tau- 



functions can be shown to satisfy the 1^-constraints (see [Q, 10|). This will be done 
in another place. 

An anti-isomorphism which can successfully be used for constructing new bis- 
pectral operators from old ones is the formal conjugation: b{x) = x, b{dx) = —dx- 
Instead of formulating a general result we shall give an example. 

Example 2.5 Let /? G C^, 7 G C^, A + 7i = - 1. Define B = {x^,Lp,Dx), 
B' = {z^ , L^, Dz), where the operators L^, Ly are Bessel operators (see or 
Example ^ (hi)). Let b be the formal conjugation: b{x'^) = z^ , b{Dx) = —D^ — 1, 
b{Lp) = {—l)^Lj. As acts locally nilpotently on B' (since it is bispectral) one 
can define bi = e~'^'^^ o b. Then we have 

biiDx) = -D,-l + L^, 

6i(L^) = (-l)^L^, 

61 (x^) = M^ = e-'^^~'^->{z^). 

My can be computed easily applying 61 to both sides of the identity 

L^x^ = {Dx + N-pi)---{Dx + N- Pn). 
In the simplest nontrivial case A = 3 we have 

= -L2 + (sDz -9-Y,li) ^7-3^'+ (9 + 2 ^ 7i) L>,- (9 + 3 7i + ^ 7,7^) +z' 
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The corresponding function ip{x,z) satisfies 

{D, + D, + = L^V = (-1)^^/3^', 

z^il, = MpTp. 

For "generic" (3 (see @) botli spectral curves are C. 

This construction can be generahzed considering bi = e^'^^^^''^ o b for arbitrary 
polynomial p. □ 

Example 2.6 The last example deals with the non- generic /? = (— 1, 1,3). In this 
case the algebra = {x'^,x^ , D^, Lfj) is larger than for generic /3 (cf. Example 
(ii))- We have the anti- isomorphism 

2 2 3 3 

X Z , X Z , 

b: D,^-D,- 1, 

1/3 1-^ 1/3. 

Let hi = e^~^/^^^'^^i^ o b. We obtain the bispectral operator 

^ ^ g{-i/3)adL^(^2) ^ ^-4p ^ -1){D- 2){D - 4) - 2x-\D + 1){D-1) + x\ 

which can be written as 

L = {d"^ - x-'^d - x){d^ + x-'^d - x). 

It coincides with a Darboux transformation of the Airy operator d"^ — x - see 
Example 5.8 for do = 1, a = oo, A = /i = (cf. also Q). 

The spectral algebra has rank 2, the spectral curve is C with a cusp at 0. □ 

We hope the above examples help in unifying the Bessel and Airy cases of which 
although similar still had important differences. 

3 Bispectral Darboux transformations 

Up to now the most widely used method for constructing solutions to the bispectral 
problem is the method of Darboux trasformations introduced in the pioneering work 
of Duistermaat and Griinbaum (see also |]l^, |l^, ^, ^ ^ ) . Here we give a general 
definition of bispectral Darboux transformations in the context of bispectal triples 
introduced in Sect. 1. 

Definition 3.1 Let {B^,B'^,b) be a bispectral triple. Let L G be a bispectal 
operator which can be factorized as 

L = QQ-^P with P,QeB^,ee JC^. (3.1) 

We call the operator 

L = PQe-^ (3.2) 
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a bispectral Darboux transformation of L and the function 

4^ = Pip (3.3) 
a bispectral Darboux transformation of ip. □ 

Theorem 3.2 The operator L is bispectral. More precisely, if f = b{L) € /C^ then 

L^ = fii (3.4) 

and 

Alp = elp, (3.5) 

w/iere 

A = 6(P)6(Q)ri. (3.6) 



Proof. Let A = ^(6*). As has no zero divisors (|3.1|) imphes 

= PL-^Q, L-^ G 1^. 
Applying the anti-isomorphism b we obtain 

A = bie) = b{Q)f-\p). 

This shows that A is a bispectral Darboux transformation of A. From ( |3.2[) and 
( ^ it follows that Ip satisfies (jsj) and (^. □ 



Theorem 3.2 is an abstract version of our bispectrality theorem from Q. How- 
ever, it formalizes only partially the arguments of this theorem. The main difference 
is that here we do not touch the subtle (in general) question when a factorization 
of the form ( p.l| ) is possible. Nevertheless, Theorem can be a part of a general 
theory of bispectral operators. Below we illustrate it by examples which have not 
appeared elsewhere. 



Example 3.3 (i) Let {B,B',b) be the bispectral triple from Example 1.4 (i), i.e. 
B = B' he two copies of the g-deformed Weyl algebra. Then L = 3^ is a bispectal 
operator with ■ip{x,z) = expg{xz), i.e. 

dg expg(xz) = expg(a;z). 

One can easily check that 

L={{a + q^x)dq - l)((a + x){a + qx)y^{{a + x)dq -q + l), 



which gives the factorization of ( |3.1| ). Then if we put 

ipi{x, z) = {{a + x)dq - q + 1) expg{xz), 

Theorem gives that 

((a -I- x)dq - g -I- l)((a -I- q^x)dq + q^ - l)((a + x){a + qx))~^'ipi{xz) = z'^ipi{xz), 
{-z{a - dq) - q + l)(-z{a - q^dq) + q^ - l)z^'^ipi{xz) = [a + x){a + qx)ipi{xz). 
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(ii) Take the same bispectral triple but put 

L = dl = x~'^{D - l){q'^D - l){q'^D - 1). 

For {a,6,c} = {0,1,2} L can be factorized as {q^^+^D -l){q-^+^D -l)x-^{q-''D - 
1). The procedure of Theorem 5]1 gives the operator 



L = x-^iq'^D - l)iq-^D - l){q-^-^D - 1). 
It is clear that repeatedly applying this procedure one can get 

L = x-^{q-'^''D - l){q~^"'~^D - l)(g^3n-2^ _ 

with any k,m,n > 0. Notice that L is a g-Bessel operator defined in Example |1.4| (ii). 
(iii) Let {B, B', b) be the bispectral triple defined in Example (ii) and put 

L = L^,, = x-2(D-/i)(L»-/2). 

Let P2 — Pi = 2a, a G Z>o. 

Then one can check that L factorizes as L = QO^^P with 

Q = g22)(l + ag2V") -g^2(l + a2:2"), 

P = (l+ax2")L»-g^ni + a9^"^^"), 
e = x2(l + aa;2")(l + ag2°2;2°), 



for a € C. Applying the recipe of Theorem |3.2| one gets easily the bispectral operators 
L, A : 

L = PQe-\ A = b{P)b{Q)f-\ 
satisfying (3^, |3.5D , where j = and 

h{P) = D{l+aL")-q^^{l+aq^"L"), 

b{Q) = (l + a(?2"L")g2D_g/32(i + aL"). (3.7) 

□ 

We end this section with the following observation. Suppose starting with a bis- 
pectral operator one produces a sequence of bispectral Darboux transformations la- 
beled by some integer parameters ni, . . . , nfc. Suppose that all the operators Lni,...,ni, 
are of the same order and that their coefficients depend rationally on rii, . . . ,71^. 
Then considering ni, . . . , as arbitrary complex parameters we get a family of bis- 
pectral operators. The standard example is from Q where from one produces 

— n{n — l)x~^ and putting n(n — 1) = c one gets arbitrary classical Bessel opera- 
tor. Another example is provided by Example |3.3| (ii). The following is a non-trivial 
example illustrating this method. 

Example 3.4 If '4){x, z) is the ordinary Airy solution of the bispectral problem, i.e. 

{dl - x)ip = zip, 

xip = {dl - z)ip, 
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then D^^ = d^id^ - x)^) = {x-^D^iD^ - 1)(D^ - 2) - - From this it is 
easy to see that ^ is the joint eigenfunction of Example \2.5\ for (3 = (0, 1,2). For 
arbitrary P € denote that function by ipp. Then we have the following Darboux 
transformation: 

V'/3i+i,/32+i,/33-2 = -(^Li3^^i32{x,drc) - a:;^V'/3l,/32,/33• 
For example the identity 

^-^^/3i+l,/32+l,/33-2(a;,9x) - Dx- l^^/3i+l,/32+l,/33-2 

is a straightforward consequence of the equation 

(-^/3i,/32 - ^)i^l3iA,l33 - D -2) = (L/3^4.i,/32+i,/33_2 - D - l)(L/3^,/32 - x). 

In this way we obtain the bispectrality of Example for /? = (i, 1 + j, 2 — z — j), 
where i,j are arbitrary integers. Because of the polynomial dependence on /3 of 
the above differential operators we can continue the joint eigenfunction for arbitrary 
/3 € such that /9i + /32 + /^s = 3. Note that in Example |2.5| one can normalize 
/3i + /32 + /33 by multiplying -0/3 by x^z-" for c G C (cf. g). □ 



^2'0/3i+l,/32+l,/33-2 



The last Example together with Examles 2^ and shows that there is a connection 
among all methods presented here. 
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